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In flat geometries, two dimensional hexatic order has only a minor eflfect on capillary waves on 
a liquid substrate and on undulation modes in lipid bilayers. However, extended bond orientational 
order alters the long wavelength spectrum of these ripples in spherical geometries. We calculate 
this frequency shift and suggest that it might be detectable in lipid bilayer vesicles, at the surface 
of liquid metals and in multielectron bubbles in liquid helium at low temperatures. Hexatic order 
also leads to a shift in the threshold for the fission instability induced in the later two systems by 
an excess of electric charge. 
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One of the main predictions of theories of disloca- 
tion and disclination mediated melting in two 
dimensions is that the transition from solid to liquid 
can be a two-stage process First, at a tempera- 

ture T = Tm, dissociation of dislocation pairs drives 
the transition from a crystal to an intervening hexatic 
phase. This liquid-crystalline phase has no translational 
order but still exhibits extended orientational correla- 
tions. Then, its quasi- long-ranged orientational order is 
destroyed in a second transition which takes place at a 
higher temperature T = Ti. Here, an unbinding of discli- 
nation pairs finally produces an isotropic liquid. 

The hexatic phase predicted by theory has now been 
observed in free standing liquid crystal films Q , in vari- 
ous types of two dimensional colloidal crystals in 
magnetic bubble arrays 1^ and in Langmuir-Blodgett 
surfactant monolayers . Despite the modest time scales 
available even on the fastest computers, there is now ev- 
idence via computer simulations for continuous melting 
and a narrow sliver of hexatic phase for hard disks |^ and 
a more substantial hexatic phase for particles interacting 
with a repulsive 1 /r^^ potential Q . 

In the above experiments, two dimensional order was 
typically probed via diffraction or by direct measure- 
ment of correlation functions in real space. It is difficult 
to use these methods when hexatic order is present in 
a curved, spherical geometry. Examples where hexatic 
order might be present include "liposomes", i.e., closed 
vesicles composed of lipid bilayers the surface of 

liquid metal droplets confined in Paul traps |pd]| , and 
multielectron bubbles submerged in liquid helium [p^ . 
Two dimensional planar multilayers of the phospholipid 
DMPC have many similarities to free standing liquid 
crystal films, which are known to have hexatic phases 
Celestini et al. have found evidence from computer simu- 
lations for extended orientational correlations at the sur- 
face of supercooled liquid metals |lj] . Electrons trapped 
on the surface of liquid helium by a submerged, positively 
charged capacitor plate have long been used to investi- 
gate two dimensional melting [|5l . Although these exper- 



iments are consistent with dislocation mediated melting, 
evidence for a hexatic phase and disclination unbinding 
remains elusive. Multielectron bubbles result when large 
numbers of electrons (10^ — 10^) at the helium interface 
subduct in response to an increase in the anode potential 
and coat the inside wall of a large sphere (10-100 micron 
radius) of helium vapor. 

One might hope that bond orientational order in a 
membrane or interface could be detected by its effects 
on the dynamics of undulation modes or capillary waves. 
Unfortunately, hexatic order couples only to the Gaus- 
sian curvature p6[ |, which vanishes for a simple sine 
wave deformation of a flat membrane or an interface 
(cf. Eq. ( [Til ) below). The situation is different, how- 
ever, when these excitations are superimposed on a non- 
trivial background geometry such as that of a sphere. 
Here, we determine the effect of hexatic order on, e.g., 
the quadrupole and octopole undulation modes and capil- 
lary wave excitations for the spherical systems described 
above. The frequency shift is large for liposomes with 
hexatic order and might also be observable in liquid metal 
droplets and in multielectron bubbles in helium. 

We first consider liquid droplets with surface hexatic 
order, either due to surface ordering or to coating with 
surface active molecules. The equilibrium shape mini- 
mizes a free energy Fd given by 

Fd = a J cIA + ^KaJ dA D,nWn, , (1) 

where a denotes the surface tension of the interface con- 
necting a droplet with density pi surrounded by vapor 
with density p^. For a general manifold with internal 
coordinates x = {xi,X2), the surface element is given 
by dA = y/^d^x, where g{x) is the determinant of the 
metric tensor gij (x) . For an undcformed sphere with ra- 
dius i?o, X = {0,^p) with polar coordinates 9 and ip and 
dA — Rq sin OdOd(p. Finally, n is a unit vector in the tan- 
gent plane with n^n* = 1 which identifies (modulo 27r/6) 
the long-range correlations in the hexatic bond directions 
]l6[ . Here, Din^ = g^^DiUk, where g^^ is the inverse of 
gij and Di denotes a covariant derivative with respect 
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to the metric gij. Thus, DiU^ = diU^ + T-'j^^n^ , where 
the are ChristofFel symbols of second kind, see e.g. 
0. The hexatic stiffness Ka ~ Ec{^t I aoY , where 
is the translational correlation length, ap is the particle 
spacing and is the dislocation core energy. The ra- 
tio KA/ksT jumps from an universal value 12 /it to zero 
when the hexatic melts into an isotropic liquid at T = 

The fundamental assumption which underlies the hex- 
atic free energy discussed above is that the configuration 
of minimal elastic energy corresponds to a vector field ng 
where no (2; -I- dx) can be obtained from no(a;) by parallel 
transport of no- On a sphere however, curvature intro- 
duces "frustration" since parallel transport of n along 
closed loops on the surface leads to a rotation of ft. For 
geometries frustrated by a nonzero integrated Gaussian 
curvature G, the state of minimum energy always has 
topological defects |T^ ]. 

These defects reduce the elastic energy by screening the 
Gaussian curvature. This point can be made more pre- 
cise by introducing a local bond-angle field 0^ the angle 
between n and some local reference frame. The singular 
part of 9 is then connected with the disclination den- 
sity [ p^p^ and the hexatic contribution F}^ to Eq. (|l|) 
becomes 119(1 



Fh ^--Ka I dA I dA' [G{x) ~ six) 



[Gix')~six')] 



(2) 



Here, 1/A is the inverse Laplacian, G{x) the Gaussian 
curvature and s{x) the disclination density |p^ . 
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with Nci disclinations of charge qi = ±27r/6 at positions 
Xi . The defects minimize F^ by arranging themselves to 
approximately match the Gaussian curvature. For low 
temperatures, large core energy Ec and spherical geome- 
tries we expect Nd — 12, corresponding to 12 5-fold 
disclinations at the vertices of an icosahedron. 

To investigate the influence of hexatic order on 
droplets, we study deformations about the equilibrium 
configuration. We expand the free energy in a small 
time-dependent displacement field SR{x,t), where R' — 
Rq + SR is the deformed surface and Rq is the radius vec- 
tor of a sphere. Here, the displacement field can be cho- 
sen to be purely normal, 6R(x,t) = Ro({x,t)N , where 
N = Rq /Rq is the normal vector of the sphere and C can 
be expanded in terms of spherical harmonics 



ax,t) = Y. E 

lm(x). 



(4) 



In the absence of defects, the expansion of F in C would 
be straightforward. On the sphere however, one has to 
deal with a distribution of discrete disclination charges 
which produces a small static icosahedral surface defor- 
mation. We initially neglect this discreteness and show 
afterwards that the corrections arising from the discrete 
nature of s{x) are irrelevant for the oscillation frequencies 
Lu{l) with I < 6. These considerations can be made more 
precise by expanding s{x) in terms of spherical harmonics 



^ 00 I 
s{x) = Go + j^Yl E 

1 = 1 m=-l 



SimYi,„(x), 



(5) 



where Go = l/i?o and sim = J^fJiU'^imi^i)- In the 
following, we first assume sim — for I > and then 
discuss the corrections due to nonzero sim- 

For Sim — with I > the hexatic order has no influ- 
ence on the droplet shape and the equilibrium configura- 
tion is a sphere with a radius Rq. However, the presence 
of hexatic order and a nonzero stiffness constant Ka does 
affect the fluctuation spectrum lu{1) of a hexatic liquid 
droplet. To see this, we solve the Navier-Stokes equation 
for the fluid in the presence of an undulating interface 
with boundary condition pO|| 
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a$(r) 



dt 
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r=(fio-l-Cflo)- 



dt 



r={Ro+CRo) + 



Ap{x), 

(6) 



where $(r, a;,t) is the velocity potential. Eq. (|^) 
thus relates the pressure difference between the in- 
side and outside of the droplet with the general- 
ized pressure discontinuity Ap which is caused by the 
shape displacement. $(r, is here given by |2^] 

'^{r,x,t) = El™ A™ Yin, (x) i^y^' for r > Ro{l + C) 

and ^r,x,t) = Ei.™ A™Yi„,(x) (^^)' for r < Ro{l + 

C). The displacement field ( and the velocity poten- 
tial <i> are related by the additional boundary condi- 
tion _RoC = RodC/dt = d^/dr\^^j^^^-^_^Q, which yields 
= -Rlrim{t)/{1 + 1) and A<^{t) = Rlh„r{t)/l. 
Upon setting Ap{x) = E;,m(^P/m(t))Yini(x), one then 
finds 



(7) 



ApUt)=\j + f^]Khm{t), 



where the generalized surface pressure of the displaced 
surface is given by Apim = —5F{ri,-a) / Ra^Tun- Since the 
interfacial energy contribution to Eq. is (cf. e.g. [pO[) 



F^ = \aRlY,\rim\\l-m + 'i) 



(8) 
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1=0 m=-l 



and the hexatic free energy reads pT[ 
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we set r;^ 
Pv < Pi) 
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^— 1 ra— — l 



\rir. 



l{l + l) 



(9) 
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Note that ijj(Z) vanishes for / = 1, corresponding to 
translations of the droplet as a whole. Eq. (^0|) also shows 
that hexatic order only affects capillary waves in a curved 
geometry: In the flat space limit of large Rq and I ^ 1 
with k = 1/ Rq fixed, one has 



fc3 
Pi 



Ka 
Rl 



(11) 



The hexatic contribution drops out as i?o ~> oo and we 
recover the result for capillary waves of a flat fluid sur- 
face Thus, it is essential to study deformations of 
a curved geometry to reveal the presence of hexatic or- 
der. In general, the undulation frequency Eq. de- 
pends on the ratio which for hexatic order at the sur- 
face of supercooled hquid metal droplets is Ka/(tRq ~ 
{^T / Rof{Ecl cra^)- This ratio becomes of order unity 
when ~ Ro- 

Nonzero coefficients sim with I > affect the mean 
curvature H{x) of the stationary droplet via the extremal 
equation |21 



2H = 2Ha 



Ka 1 



E 



(/-!)(/ + 2) 



aRlRo 1(1 + 1) 

^' l — l rn— — l 



Slm^h 



with Hq =^ I/Rq, leading to static surface deformation 
coefficients = simKA/o'Rol{^'^^) which vanish in the 
limit Rq ^ oo. Thus, for Ka 7^ the defects effectively 
repel each other and deform the droplet, cf. Fig. |^. How- 
ever, as we show below, nonzero sim have no influence on 
the frequencies uj{l) for < ^ < 6. 




FIG. 1. Shape of a liquid droplet with hexatic order and 12 
disclinations lying on the vertices of an icosahedron (surface 
deformations associated with the defects are exaggerated) . 



Next, we discuss the influence of hexatic order on the 
fluctuations of a vesicle. Here, the free energy is given by 
00 



Fy^^K J dA {2Hf + '^G y dAG 



Ka / dA D,nW'n-i, 



(12) 



where k and kg are the mean and Gaussian rigidity, re- 
spectively. In the following analysis, the second term of 
Eq. ( |l^ ) can be neglected since we only consider surface 
shapes which are topologically equivalent to a sphere. 

Since k plays a similar role for vesicles as a plays for 
droplets, hexatic order should lead here to similar effects 
on the spectrum io{l). However, for vesicles the fluctua- 
tions are overdamped and one has to analyze the Stokes 
equation. This can be done by generalizing the approach 
of 1^ to hexatic membranes. One then obtains |2^ (ne- 
glecting for simplicity the volume constraint) 



'^(0 = -*^(/-l)(^ + 2) 



kI{1 + 1) + Ka 



l{l + l) 



(13) 



where 77 is the liquid viscosity and T{1) = l{l + 1)/(2Z -|- 
1)(2P + 2Z — 1). In the flat space limit one has 



Kk + o k 



(14) 



in agreement with ||23| for i?o ^ 00. 

The frequency shift (^3|) now depends on the ratio 
Ka/ n. However, as i?o — ^ 00, we expect that Ka — 4k 
(a universal result for flat hexatic membranes at long 
wavelength |2^) leading to a frequency enhancement by 
a factor ~ 13/9 ~ 1.44 for the I = 2 quadrupole mode. 
Thus, bond orientation order has a strong effect on the 
fluctuations of a membrane and should have experimen- 
tally observable consequences in dynamical light scatter- 
ing (2|. 

For vesicles, the presence of a finite number of de- 
fects also leads to an equilibrium configuration with a 
deformed surface. The mean curvature H(x) of the sta- 
tionary vesicle now satisfies 



Ka 1 



2i7 = 2ifo + ^^E E 



K Rq 



l—l m- 



(^-l)(^ + 2) 
P{1 + 1)^ 
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leading to nonzero coefficients 



si^Ka/kI\1 + 1)2 



in the ground state. However, icosahedral symmetry in- 
sures that Sim = unless / = 6, 10, 12, so corrections 
of order sim have no infiuence on the frequencies lo{1) 
for small I. Provided the positions of the discfinations 
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remain fixed in the lipid matrix on the time scale of an 
undulation {uj{l = 2) ~ 40Hz for a Ifi vesicle) the dis- 
persion relation ( p^ remains valid for < ^ < 6. Be- 
cause disclination motion is catalyzed by absorption and 
emission of dislocations with spacing ^t, the disclina- 
tion diffusion constant is « {ao/^T)^Diipid, where 
Diipid ~ 10~^cm^/sec. The estimate {aol£,TY — 10^^ 
suggests only minor disclination motion during an undu- 
lation period. Similar estimates show that disclination 
motion is negligible during a capillary wave period for 
hexatic droplets. 

Finally, we will discuss niultielectron bubbles in liquid 
''He. These bubbles can undergo both a freezing tran- 
sition and a shape instability. Thus, here hexatic order 
affects not only the fluctuation spectrum but also the 
instability-threshold for fission. The free energy of a mul- 
tielectron bubble — Fd + Fc is that of a droplet (cf. 
Eq. (|l|)) with an additional Coulomb contribution 



dA 



I P{x)p{x') 



(15) 



where p{x) denotes the charge distribution on the surface 
and e is the '*He dielectric constant. In an equilibrium 
fluid, p — eN/AirRQ for a sphere with N electrons. 

Within the approximations described above, one now 
flnds (with I > and neglecting the density inside the 
bubble) ll 



X {I 



{i-m + i) 



2) -4. 



, J- i-^^ 



Ka (/-l)(/ + 2)2 



cri?5 



(16) 



where R^^. = (eN)'^ /IGnae is the critical radius for mul- 
tielectron bubbles without hexatic order |25|. Thus, for 
Ka = spherical bubbles become unstable to flssion if 
R < Rcr, i.e. uj'^{l = Ic) <0 for R < R^. and Ic = 2. For 
Ka 7^ the stability of charged bubbles is enhanced by 
the hexatic order of the electrons on the sphere. Thus, 
for Tm <T <Ti one stiU has Ic ^ 2 but u'^Hc = 2) < 
for Ro < R'c with R'^ < Rc. The icosahedral symmetry of 
the deformed shape with s;,„ ^ is too high to have an 
influence on the flssion instability which occurs at Z = 2. 

Because the electrons (which determine Ka) are far 
apart relative to the helium atoms (which determine 
a) Ka/<jR^ will be smaller than for droplets of super- 
cooled liquid metals with Rq ~ Rc. For helium-bubbles 
with N ~ 10^ one has Rc — 10/im and we expect that 
~ 10"'^. Charged metal droplets undergo the same 
flssion instability. Thus, it might be possible to detect 
hexatic order by investigating the stability of charged 
liquid droplets in Paul traps. 

The effects discussed here are even larger if the hex- 
atic phase is bypassed and one freezes directly into a 
curved two-dimensional solid with shear modulus /i pw . 



The resulting frequency shifts can be estimated by set- 
ting Ka — P-Rq in the formulas above (although the l- 
dependence will be different). There are also interesting 
consequences of hexatic order for dynamics of liquids and 
membranes in a cylindrical geometry pl[ | . 
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